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A METHOD FOR FINDING THE ROOTS OF THE EQUATION 
/(x)=0 WHERE / IS ANALYTIC 
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ABSTRACT 

An expression is derived giving the roots of equations involving only functions 

which may be expressed in power series. The expression takes the form x n =x — 

11 JSL / (r) (x ) 

H*zi, w here T^ (— 1) **«' &»-<=* and a r = J ., ° > and x n approaches the value of 

the root as n increases without bound. Methods are given for the treatment of 
certain special forms. An example is given of the application of the method to 
the calculation of real roots. 
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I. INTRODUCTION 

The use of recursion formulas in evaluating roots of equations was 
described by Euler. 1 Since that time the formulas and the related 
determinants have been presented by various authors, 2 the most 
recent example being the series involving determinants developed by 
E. T. Whittaker. 3 However, three of the quantities necessary to the 
evaluation of the r th term of this series are sufficient to determine the 
sum of the first r terms, so that the series form is not the most conveni- 
ent one for calculation. In the present paper the recursion formulas 
are developed in an elementary way, and put in a form which is con- 
venient for calculation. 

1 Introductio in Analysii Infinitorum, 1, chap. 17 (1748), or the translation by J. B. Labey, p. 257 (Paris, 
1796). 

2 F. Cajori, A History of Mathematics, p. 227. (McMillan, London, 1919.) 

3 Proc. Edinburgh Math. Soc. 36, 103 (1918), or Whittaker and Robinson, Calculus of Observations, 
p. 120. (Blackie and Son, London, 1924.) 
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II. DERIVATION OF EQUATIONS 

Eequired the value of x, which is a root of the equation f(x) = 0, 
where f(x) is analytic in the neighborhood of x. Then if an arbi- 
trary approximation to x, x , lies in the neighborhood of x the function 
f(x) is given by the series 

f( x )=a +a 1 (Ax)+&(Axy+ .... +%(Axy+ (1) 

where the coefficients a r =f ir) (x ) are given and Ax=x—x . 
Define 

g n{x) = l +Kl (& x ) + b -^( Ax y- + + A^i_ (A3 .)»-l (2) 

where the n— 1 constants bn.i,b n ,2, ,&*,»-i are so chosen that n-1 

consecutive derivatives (beginning with the second) of the function 
<p n (x) vanish at x , where 

<P*(x) =f(x)g n (x)==<p n (x ) + (Ax)<p n '(xo) + ( ?l ^ 1 )i ^ (w+1) fa) + 

Then 

<Pn(xo) = 0=a 2 + 2a 1 b n>1 + a b n , 2 

<p"n (%o) = = ^3 + 3a 2 6 w ,i + 3ai6 ni 2 + a Q b n , 3 



^nM = = a T +\lja r ^ibn t i+ +( r r _ l )aib nj _ l + a b nj ^ 



^l\x )^0=a n +(^)a n _ l b n ,i+- • +L^2)^ 6n,n - 2 ^(n— l) ai6n,n - 1 

The symbol ( . ) is used to represent the number of combinations of 

r things taken i at a time. 

The equations 4 may be solved for <Pn(x ) y with the result 

*«*)"-§; (5) 
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where A n is the cofactor of a x in the determinant of the n th order 
having for the element of the r th row and c th column ( - la r _ c+ i. 
The determinant may be written 



D n = 



as 



a 
2oi 
3a,2 

4a 3 



. . 


. . . . 


. . 





do . . 
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. . 





3ai . . 
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. . 





6d2 


. . . . 


. . 






(i)*- 1 G) ar - 2 (c-iK 



/n-l\ fn-i\ /^-1\ /w-l\ 

On-i (^ j Ja n _ 2 (^ 2 )a n - 3 . {<._! )an-c . . ^ n _ 2 / 11 a ° 

(1)0,-1 (2)^-2. ( c * X )<Wi • • ( n ^ 2 ) a2 (n-l) ai 



Equation 3 may be put in the form 

<p n (x) =a + (Ax) <p n ' (xq) +R n 



(6) 



where R n represents the sum of the terms containing powers of 
Ax higher than the n th . Solving equation 6 for Ax and putting 
<Pn(x)=0 } we have 



Ax= 



a +Rn 



<Pn'(Xo) 



We may write 



A n x=- 



do 

<Pn(Xo) 



(7) 



(8) 



Then, as lim R n = 0> lim A n x=Ax. 

n— > 00 n— > 00 

From equations 8 and 5 we have 



A„x= 



a A n 



D n 



(9) 



It may be shown that ^L w =ni? ft _i. Multiplying each element of 

r+1 /r+l\ 

Dn-i by 9 the element becomes f ' J a r _ c+1 , which is the element 

of -4 n . In this process D n ^i is multiplied by , ZTT\Y or n " Thus 

A n =nD n ^ (10) 



598 Journal of Research oj the National Bureau oj Standards [ vol. u 
Then from equations 9 and 10 

AnX= -^=i (11) 

A recursion formula for D n may be obtained by expanding it on 
the elements of the n th row. The result is 

^(-lyQy&O^Dn^O (12) 

This ma} r be simplified by writing U r =-rr and a r ss-~ 

r la o r \a 

Equation 12 then becomes 

S(-l)W»-i=0 (n=l,2,3, . . . .) (13) 

i=0 

Now equation 11 may be written 

A„x=-%f (14) 

Clearly U =D =1, in this development, but in general the value of 
U may be chosen arbitrarily and another value may be more con- 
venient. 

The product of k consecutive convergents to Ax is 

r=n-\-k tj 

n Ar x=(-l)* 7 p- (15) 

r=n+l Un+k 

The limit of this product as n increases without bound is (Ax) k . 

III. TREATMENT OF SPECIAL FORMS 

In c&sef(x) = F(x 2 ), and x is taken equal to zero, the derivatives of 
odd order are zero, and the functions U 2m +i are also zero. Thus, 
A2 m x=0 and |A 2w+ ix| = ». Taking n even and k=2 in equation 15 
gives 

(A 2m+1 x)(A 2m+2 x)=-ff^ (16) 

This expression may be evaluated, giving a convergent to (Aa;) 2 . 
For instance, let / (x) = cos x=Q, and take x =0. Then U 2m becomes 

W 
tt^v)' where E m is Euler's number of order m, or Bernoulli's number 

of order 2m, and 

(^(^#)= (2m+1 f m+2) ?~ (17) 
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Now lim(A 2m+1 x) (A 2m+2 x) = (Ax) 2 and fr J2m+l)(2m+2)# w = ^ 

so that (Ax) 2 =-j or Ax= ±«* which gives the roots ±^ for cos x=0. 

The case in which j{x) = F[(x—~c) v ] and x Q =c may be treated in a 
different way. Put (x— c) p = z, and find Az for z Q =0. Then Ax= 
(As) 1/2 \ 

IV. APPLICATION OF RESULTS 

The familiar equation x 3 — 2x— 5 = does not afford a good test of 
the method of this paper. If x is taken as 2, the values of a i9 <x 2y 
and a 3 are —10, —6, and —1. Thus, the successive values of U n are 
integers, and may be written down immediately as far as U 4 , which 
is sufficient for the calculation of the root correct to six decimal 
places. 

The advantages of the method are better illustrated by the solu- 
tion of a transcendental equation, such as 



/(x)=x+sin x+log x + e z — 5 = 
Taking #0=1.1, we find 



(18) 



a = 0.090 683 6 
ai = 5.366 853 
a 2 = 1.286 51 
a z = 4.053 
a 4 =-0.2 


ao/2 = 0.643 26 
a 3 /6 = 0.675 5 
a 4 /24=-0.01 


Table 1. — Table of the products m 


cessary for the calculation of 





ao=l 


a\ = — 


«2 

2«o 




a4= 24^ 


Ui 

u 2 

Us 


1 
59. 182 18 
3, 495. 437 
206, 455. 2 
12, 194, lie 


59. 182 18 
3, 502. 530 
206,867.6 
12, 218, 46 9 


7. 093 4 
419. 80 
24, 794 


7.449 
440.8 


-0.1 























J7i=a! A 1 x=- 1/C7i= -0.016 90 

TJ 2 =a x U x -a 2 A 2 x=-£/ 1 /£/ 2 = -0.016 931 3 

U z =a l U 2 -a 2 U l + a s A 3 x=- U 2 /U z = -0.016 930 7 

Di=ai'Di— otJJt+aiUi—cn A 4 x=-£7 3 /£/ 4 = -0.016 930 7 

It may be seen that the value of a 4 was not significant in the calcula- 
tion of Z7 4 . The accuracy of the calculation is not increased by the 
last step, as A 3 x contains no error which is significant with respect to 
the value of a used. Then x 3 =x 4 = 1.083 069 3. The value of the 
root correct to eight decimal places is 1.083 069 29. 

It should be noted that in solving equation 6 for Ax, <p' n (x ) was 
used as a divisor. Thus, if <p' n (x ) is zero or infinity, as may be the 
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case for certain values of x and n, equation 7 does not hold. Gen- 
erally this difficulty may be avoided by taking a higher value of n. 
<p' n (xq) will not ordinarily be zero, however, if x is a good approxi- 
mation to a root. 

As the derivation of equations 13 and 14 does not involve any 
assumptions other than that/(x) is analytic, the equations hold for 
complex values of the coefficients or the roots. In the evaluation of 
complex roots of real equations it is convenient to use a complex value 
for Xq, as equation 15 must be used if x is real. This treatment is not 
satisfactory, as the sequence U n /U n + k does not usually converge 
rapidly in such a case, and also because extraneous roots may be 
introduced. 

V. SUMMARY 

An expression is derived giving the roots of equations which involve 
only functions which may be expressed as power series. The expression 
involves only the values of the function and its derivatives at a given 
value of the variable. 

Methods are given for the treatment of certain special forms. 

An example is given of the application of the method to the evalu- 
ation of real roots. 

Washington, November 22, 1934. 



